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Abstract. Dark energy with the usually used equation of state p — wp, where 
w — const < is hydrodynamically unstable. To overcome this drawback we consider 
the cosmology of a perfect fluid with a linear equation of state of a more general form 
p = a{p—po), where the constants a and po are free parameters. This non-homogeneous 
linear equation of state provides the description of both hydrodynamically stable 
(a > 0) and unstable {a < 0) fluids. In particular, the considered cosmological 
model describes the hydrodynamically stable dark (and phantom) energy. The possible 
types of cosmological scenarios in this model are determined and classified in terms 
of attractors and unstable points by the using of phase trajectories analysis. For the 
dark energy case there are possible some distinctive types of cosmological scenarios: 
(i) the universe with the de Sitter attractor at late times, (ii) the bouncing universe, 
(iii) the universe with the Big Rip and with the anti-Big Rip. In the framework of a 
linear equation of state the universe filled with an phantom energy, w < — 1, may have 
either the de Sitter attractor or the Big Rip. 



§ dokuchaev@inr.npd.ac.ru 



Dark energy cosmology with generalized linear equation of state 



2 



1. Introduction 

The astronomical observations indicate that the expansion of our Universe accelerates 
[1]. In the framework of the General Relativity this means that about two thirds of 
the total energy density of the Universe consists of dark energy: the still unknown 
component with a relativistic negative pressure p < — p/3. The simplest candidate for 
dark energy is the cosmological A-term or vacuum energy. During the cosmological 
evolution the A-term component has the constant (Lorentz invariant) energy density p 
and pressure p — —p. However the A-term requires that the vacuum energy density be 
fine tuned to have the observed very low value. For this reason the different forms of 
dynamically changing dark energy with an effective equation of state w = p/ p < —1/3 
were proposed instead of the constant vacuum energy density. As a particular example 
of dark energy, the scalar field with some a slow rolling potential (quintessence) [2] 
is often considered. The possible generalization of quintessence is a /c-essence [3], the 
scalar field with a non-canonical Lagrangian. 

Present observation data constrain the range of equation-of-state of a dark energy 
as —1.38 < w < —0.82 [4]. Recently the Chandra observations [5] gave the similar 
restrictions on the value of w. These data do not exclude the possibility of our universe 
to be filled with a phantom energy [6]: the energy with a super-negative equation of 
state w < —1 (note, however, that the dark energy with the equation of state w evolving 
from the quintessence-like w > — 1 in past, to phantom- like w < — 1 at present, provides 
the best fit for supernova data [7]). The different aspects of phantom cosmology were 
considered in [8]. The phantom energy is usually associated with the phantom or ghost 
fields — the scalar fields with wrong-sign kinetic term. It is known that such fields are 
unstable due to quantum instability of vacuum, unless the kinetic term stabilized by 
the ultra-violet cutoff [9]. The intriguing possibility of constructing of effective w < —1 
from the scalar field with correct-sign kinetic term was proposed in [10]. In [11] the 
authors constructed the theory with the effective -u; < — 1 in the braneworld models. 

The presence of phantom energy in the universe leads to the interesting physical 
phenomena: the possibility of the Big Rip scenario [6] , the black hole mass decreasing 
by phantom energy accretion [13] and a new type of wormhole evolution [14]. The 
alternatives to a scalar field model are the perfect fiuid models such as a Chaplygin gas 
[15]. Hao and Li [16] demonstrated that w — —1 state is an attractor for the Chaplygin 
gas and equation of state of this gas could approach to this attractor from the either 
w < —1 OT w > —1 sides. 

In this paper we analyze the perfect fluid model with a general linear equation 
of state p = a{p — po) where a and po are constants. This is a generalization of a 
homogeneous linear equation of state (po = 0) and is suitable for the modelling either 
the linear gas with p > or the dark energy with p < 0. The advantage of the use of 
the general linear equation of state is the possibility to describe the dark energy with a 
positive squared sound speed (for the usually considered equation of state p — wp the 
dark energy is hydrodynamically unstable, because dp/dp — w < 0). The considered 
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linear model is reduced to the perfect fluid model with w — const at the particular case 
of po = 0. 

In the framework of this model it is possible to find analytical cosmological solutions 
for the arbitrary values of parameters a and po- We show below that this linear model 
(unlike the Chaplygin gas model) describes distinctively different types of cosmological 
scenarios: the Big Bang, Big Crunch, Big Rip, anti-Big Rip, solutions with de the Sitter 
attractor, bouncing solutions, and their various combinations. 

The paper is organized as follows. The principal part of the paper is Sec. 2 in which 
the basic properties and particular analytical solutions for a perfect fluid cosmology with 
a linear equation of state arc derived. We show that dark energy with the equation of 
state p — ol{p — po) may be effectively reduced to the A-term and to the simplified 
linear equation of state pcs = WesPeS- In this case, either the effective density pcs or the 
density of an effective A-term may in general have the negative values. In dependence 
of signs of parameters a and po the four different cosmological scenarios are considered. 
To study the cosmological dynamics of the universe for different scenarios the phase 
plane analysis is used. In Sec. 3 we discuss the restrictions which can be set on this 
cosmological model from the astronomical observations. The discussion of the results 
and their further possible generalization is presented in Sec. 4. 



2. Lineeir equation of state 

We consider the flat Friedman-Robertson- Walker universe filled with a perfect fiuid. 
For the sake of simplicity we will call below this fluid as dark energy. Using the unit 
conventions SttG/S — c — 1 the corresponding Einstein equations for this cosmology 
can be written as follows: 

= PDE, (1) 

Pde = -3i?(pDE + Pde), (2) 

where H — d/a is the Hubble parameter and pde and pde are the energy density and 
pressure of a dark energy correspondingly. In this paper we consider a perfect fluid with 
the linear equation of state of the following general form: 

Pde = Q:(pde - Po), (3) 

where a and po are constants. Certain features of cosmology with po 7^ and < a < 1 
was considered in [12]. The equation of state (3) was used in [13] to describe the dark 
energy accretion onto the black hole. Note the difference of this equation of state (3) 
from the commonly used one p = wp with w = const. Our simple generalization allows 
to include dark (and also phantom) energy with a positive squared sound speed of linear 
perturbations cj. = dp/ dp = a > 0. 

One can use the observational restrictions [4] on the allowable range of equation of 
state, —1.38 < w < —0.82, to put the restrictions on the parameters of a linear model 
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Figure 1. Observational restrictions [4] on the allowable of the equation of state, 
— 1.38 < w < —0.82 recalculated and presented for parameters a and po of the linear 
model of dark energy. Here a and po are the parameters in (3) and Pde is the present 
density of dark energy. 



(3). Substituting (3) into the equation of state w = Pde/pde we obtain the range of 
allowable parameters parameters a and po, shown in the Fig. (1). 

The equation of state (3) can be reduced to the "effective" cosmological constant 
and the dynamically evolving dark energy by redefining the fluid density and pressure 
in the following way: 

Pde = Pa + Pa, Pde = Pa + Pa, (4) 



where 



and 



Pa = -Pa, Pa = aPa, (5) 



apo apo /«x 

PA = T— — , Pa ^ Pde-— — • (6) 

1 + a 1 + a 

In Eqs. (4,5,6) we denote Pa and p/^ as the density and the pressure of the "effective" 
cosmological constant. Correspondingly, pa and Pa are the density and pressure of the 
dynamically evolving part of dark energy. The values pa and p\ obviously obey the 
relation p\ — —3H{p\ +pa) and therefore Pa and Pa obey the equation: 

Pa = -SH{pa+ Pa)- (7) 

The cosmology with A-term and a linear equation of state (for p/p > —1) was considered 
previously in [17]. Note, however, that in our consideration either pA or pa may be 
negative (nevertheless the sum pA + Pa must be positive) and also equation of state for 
dynamically changing part of dark energy may be super-negative. This leads to the new 
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distinctive types of cosmological scenarios considered below. The signs of pa and p\ 
are conserved during the universe evolution as can be seen from (7) and (2). Below we 
obtain the full analytical solutions of evolution of the universe. In both cases p + p > 
and p + p < we find from (7) and (2): 

A 

where constant A may be either positive or negative. Two different asymptotic regimes 
are possible: for \pa\ <S |pa| the universe behaves like the de Sitter universe with 

p ^ Pa, P ~ -Pa, (9) 
and for \pa\ ^ |pa| the universe is filled with dynamical part of dark energy: 

p Pa- (10) 
Using (1) and (8) we obtain the relation between the differentials da and dt: 

±dt. (11) 



da 



a^J PK + Aa-^^^+'^) 

In dependence on signs of pa and A there can be three different results for integration 



of relation (11). For pa > and A > we find: 

-1 ~| 2/re 



a{t) 



j3 sinh 



(12) 



Here we denote k = 3(1 + «) and f3 = A/ p^\. The choice of the upper or lower sign in 
(12) depends on the sign "+" or "— " in (11). The expression (12) was first obtained in 
[17]. The asymptotic behavior of (12) for t ^ is given by 

a=(±fv^t] , (13) 
The Eq. (12) for t — > ±oo can be rewritten as follows: 



exp {±^/pX t) , 



(14) 



The expression (14) is valid only if k > 0. For k < the signs in exponent in this 
expression should be changed to opposite ones. 

In the case pa > and A < the integration of (11) gives: 

2/k 

(15) 



^(^) — { cosh 



In slightly different form the solution (15) was obtained in [18] for the restricted case 
— 1 < a < 0. The asymptotic behavior of a scale factor for t — > is given by: 



a{t) = {-/3) 



1 + 4 PA 



(16) 



For t ±oo the scale factor evolution is described by (14). Finally, for pA < and 
A > one can find the evolution of the scale factor of the universe: 



a{t) 



P sin 



2/k 



[17) 
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In the limitation — 1 < a < 1 the similar solution was found in [19] for the cosmology 
with negative anti-de Sitter A-term plus a scalar quintessential field with special form 
of a potential. For t ^ the asymptotic behavior of scale factor is given by (13). The 
behavior of a scale factor near t ~ tt/ {npiCj can be described as follows: 

2/fv 

(18) 

Correspondingly, the asymptotic behavior of scale factor at t — > ±27r / {kp\) is given by 

a= U^-p - -^/At] . (19) 

To find the attractors and unstable points of the solution of the Friedman equations we 
use the method of the phase trajectories. Denoting y = Pa +Pa we obtain the following 
system of equations 

H = -3y/2, y = -3{a + l)Hy. (20) 

All phase trajectories of the system lie on the curve — which can be rewritten 
as y = {a + l)i7^ — apo- The particular form of this curve and direction of the phase 
trajectory depends on the signs of 1 + a and apo. The four cases are possible: 
(i) 1 + a > and apo > 0. 

1 /2 

There are two singular points at y = axis. At the point {H = pj( , y = 0) the 
linearized system of evolution equations is 

H = -3y/2, y = -3{l + a)p]l'y, (21) 

The first eigenvalue and the first eigenvector of the system equal to zero. The 
second eigenvalue equals to A2 = 3(1 + a)p]l'^ < and the second eigenvector is 

1 /2 

{3/2,3(l + a)p;'^}. We can see that the considered point is the attractor. The same 
linearization method is used for the universe evolution analysis near other singular 
points. The second singular point (H — —pY'^, U — ^) corresponds to the unstable 
equilibrium state. The another interesting singular point is {H — 0, y — —apo), in 
which the universe reaches the zero density pde = and its contraction is changed to 
the expansion. The phase trajectories near the singular points are plotted in the left 
panel of Fig. 2 with the directions of evolution marked by arrows. In the auxiliary right 
panel of Fig. 2 the evolution of pressure p as a function of the dark energy density p is 
shown. 

The right branch of parabola in the left panel of Fig. 2 {y > 0, H > 0) corresponds 
to the solution (12) with an upper sign at i > 0. In this case the universe is filled 
with non-phantom energy. The universe starts from the Big Bang, corresponding to the 
initial values of density Pde,? = +00 and scale factor Ui — 0. The asymptotic behavior 
of a scale factor at a is described by (13) with an upper sign. The relation (13) 
with an upper sign describes also the evolution of a scale factor for Pa = for all t. At 
the late times the universe approaches to the de Sitter regime (9). The sign of p + p 
is conserved during the evolution of the universe. The asymptotic behavior of a for 
t — > cxD is given by (14) with an upper sign. The left branch of parabola in left panel of 
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Figure 2. Phase trajectory of the universe and diagram of evolution in coordinates 
(p, p) are shown in the case of dark energy parameters 1 + a > and apo > 0- In 
the left panel the phase space diagram of evolution in coordinates H and y = Pa+Pa 
is presented. The direction of evolution is shown by arrows. There are two unstable 

1 /2 

singular points, {H = —pfl , y — 0) and {H = 0, y = —apo), and one attractor 

1 /2 

{H = pfl , y = 0). In the right panel the evolution of is shown in coordinates p and 
p. The evolution from p = oo to p = pA corresponds to the right branch of parabola 
{y > 0, H > 0) of the phase space diagram. The reverse process (from p = p\ to 
p = oo) corresponds to the left branch of parabola {y > 0, H < 0) of phase space 
diagram. The evolution from p = pA to bounce p = and back to p = pA corresponds 
to the middle branch of parabola {y < 0) of phase space diagram. 

Fig. 2 (y > 0, if < 0) corresponds to the reverse process with respect to the examined 
above and is described by (12) with a lower sign. Note that t is taken to be negative 
in this case. The asymptotic behavior for t — > and {t — > oo) is given by (13) and 
(14) correspondingly with lower sighs in both cases. A middle part of the parabola in 
left panel of Fig. 2 {y < 0) corresponds to the solution (15). This particular case was 
considered in [18] as a simplest example of phantom cosmology without a Big Rip. The 
universe in this case is filled with a phantom energy. Because of the specifics of a linear 
equation of state, the considered universe with a phantom energy is not born in the 
Big Bang. Instead of this, the universe starts in this case from the initial scale factor 
Oj = +00. For i — > — oo it behaves like the de Sitter universe (however reversed in time) 
according to (14) with an upper sign and bounces at the minimal value of a scale factor 
flmin = i^Py^'^- In this moment the state of dark energy is very special: the pressure is 
finite but the total energy density is zero, Pa + Pa — 0. Near the bounce this universe 
can be described by (16). After the bouncing the universe expands and at the late times 
it approaches to the de Sitter state (14) with an upper sign, 
(ii) 1 + a > and apo < 0. 

The density of the A-term is negative, pa < 0. The linearized system of the 
evolution equations at the point {H = 0, y = —apo) is 

H = 3apo/2 -3{y + apo) /2, y = 3{1 + a)apoH. (22) 

Universe expands starting from the Big Bang, reaches the zero density at the point 
{H = 0, y = —apo) and then its expansion changes to the contraction. The phase 
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Figure 3. Phase trajectory of the universe and diagram of evolution in coordinates 
(p, p) are shown in the case of dark energy parameters 1 + a > and apo < 0. In the 
left panel the phase space diagram of evolution in coordinates H and y = p^ + Pa \s 
presented. The direction of evolution is shown by arrows. There is only one singular 
point, {H = 0, y = —apa), which is unstable. In the right panel the evolution is shown 
in coordinates p and p. The universe evolves from p = oo to p = 0, which corresponds 
to the right branch of parabola {H > 0) of the phase space diagram. At the point 
p = (corresponding to the point {H = 0, y = —apo) of the phase space diagram) 
there is a bounce. Then the universe evolves back to p = oo which corresponds to the 
left branch of parabola {H < 0) of the phase space diagram. 



trajectory and diagram of evolution p—p are shown in the Fig. 3 (the similar cosmological 
behavior was obtained in [19] for the universe with negative A-term and quintessence 
with special potential) . The universe is filled with non-phantom energy and starts from 
the Big Bang. The evolution of a scale factor is described by (17) and at early times 
by the asymptotic (13) with an upper sign. This asymptotic is also valid for Pa = 
at all t. During the finite time At = n/{K,p\) a scale factor a reaches the maximum 
value Umax = {—PY^'^ at which the universe bounces. Near the bounce the behavior of 
a scale factor can be described by (18). After the bounce the universe contracts and in 
a finite time At it collapses to the Big Crunch. The scale factor of the universe near 
the collapse is given by (19). 
(iii) 1 + a < and apo < 0. 

At the zero density in the point {H — 0, y — —apo) the universe expansion is 
changed to the contraction. The singular point {H — —pj^-, y = 0) is an attractor and 
respectively {H — p^l"^ , y = 0) is an unstable equilibrium point. The phase trajectory 
and diagram of evolution p ~ p are shown in the Fig. 4. 

First we consider the right branch of the parabola {y < 0, H > 0). In this case the 
universe is filled with phantom energy, p+p < 0. The cvohition of a scale factor is given 
by (12) with an upper sign. Note that t is taken to be negative in this case. Starting 
from the scale factor = the universe expands during an infinite time to the final Big 
Rip. At i — > — oo a scale factor of the universe is given by (14) with an upper sign. Near 
the Big Rip a scale factor a is described by (13) with an upper sign (note that k and t 
are both negative). The left branch of the parabola (y < 0, if < 0) corresponds to the 
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Figure 4. Phase trajectory of the universe and diagram of evolution in coordinates 
(p, p) are shown in the case of dark energy parameters 1 + a < and apo < 0- In 
the left panel the phase space diagram of evolution in coordinates H and y = Pa+Pa 
is presented. The direction of evolution is shown by arrows. There are two unstable 

1 /2 

singular points, {H = pfl , y ~ Q) and (H = 0, y = —apo), and one attractor 

1 /2 

{H = —P[l ,y = 0). In the right panel evolution of the universe is shown in coordinates 
p and p. The evolution from p = p^ to p = Poo corresponds to the right branch of 
parabola {y <{), H > Q) of the phase spac;c diagram. The reverse process from p = oo 
to p = PA corresponds to the left branch of parabola {y < 0, H < 0) of phase space 
diagram. The evolution from the p = pA to bounce at p = and back to p = pA 
corresponds to the middle branch of parabola {y > 0) of phase space diagram. 



reverse process with respect to the examined above and is described by (12) with a lower 
sign. Time t is positive during the whole evolution. In this case the universe is filled 
with phantom energy and it starts from the "anti-Big Rip" solution (with the infinite 
value of the initial scale factor), contracts during an infinite time to the final af = 0. 
The asymptotic behavior for t — > and t +oo is given by (13) with the lower sign 
and (14) with the lower sign correspondingly. The middle part of the parabola {y > 0) 
corresponds to the solution (15). The universe is filled with a non-phantom energy and 
starts from = 0. At t — >■ —oo the scale factor of the universe is given by (14) with the 
upper sign. After the bouncing at the maximum value of a scale factor amax = {—PY^'^ 
the universe begins to contract. Near the bouncing point the scale factor is given by 
(16). At i — > +00 the universe behaves like (14) with the lower sign, 
(iv) 1 + a < and apo > 0. 

In this case the density of the A-term is negative: p\ < 0. The universe reaches the 
zero density at the point {H = 0, y = —apo), where the contraction is changed to the 
expansion. The phase trajectory and diagram of evolution p — p is shown in the Fig. 5. 
The universe is filled with a phantom energy and the solution for a scale factor is given 
by (17). The universe is born in the "anti-Big Rip" state and the scale factor is given by 
(13) with the lower sign near t 0. Then the universe contracts to the minimum state 
Qmin = {—Py^'^ in a finite time At — tt/{kp\), and after bouncing it begins to expand. 
In time At the universe comes to the Big Rip where the scale factor of the universe is 
given by (19). 
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Figure 5. Phase trajectory of the universe and diagram of evolution in coordinates 
(p, p) are shown in the case of dark energy parameters 1 + a < and apo > 0. In the 
left panel the phase space diagram of evolution in coordinates H and y = pa + Pa is 
presented. The direction of evolution is shown by arrows. There is only one singular 
point, [H — 0, y = —apo), which is unstable. In the right panel the evolution of the 
universe is shown in coordinates p and p. The universe evolves from "anti-Big Rip" 
state {p = oo) to the bouncing point {p = 0), which corresponds to the left branch of 
parabola {H < 0) of the phase space diagram. Then the universe evolves to the Big 
Rip (p = oo) which corresponds to the right branch of parabola {H > 0) of the phase 
space diagram. 

Different cosmoiogicai scenarios discussed above in tliis Section are summarized in 
Table 1. 

3. Universe with dark energy and matter 

If wc add tlie usual matter (dark matter, baryons and radiation) to the considered 
universe then the differential equation for a scale factor evolution would have the 
following form: 

g)' = //o'Ko(l + Z)^ + Q^,o(l + Zf + fiA,o + Jla,o(l + zf^'^^ (23) 

where Hq is the current value of the Hubble constant, 2; = ao/a — lisa redshift, fi^o 
and flmfi — 0.3 are the cosmoiogicai density parameters of radiation and non-relativistic 
matter correspondingly. The value of the cosmoiogicai density parameter of "A-term" 
is: Qa,o — [(y/{a+ l)]po/Pc,0) where Pc,o is the current critical density. The contribution 
of the dynamically changing part of dark energy is the following: D,a,o — ^q,o — ^a,o, 
where Q^q ^ 0.7 is the quintessence density parameter. Wc suppose that the universe 
is flat with Qr,o + ^m,o + ^g,o = 1. In the case of pure A-term (a = — 1, po = 0), one 
should omit the last term in the brackets in (23) and take = ^q^- 

The observational restrictions on the w (—1.38 < w < —0.82), derived from the 
SN data, limit only the behavior of the universe at small redshifts 2; ~ — 1. In 
the case of Po/Pde < 1 (see Fig. 1) the dark energy does not change drastically the 
evolution of the universe at high redshift in comparison with a pure A-term case. The 
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Table 1. Possible cosmological scenarios in the case of a generalized linear equation 
of state p = a{p — pq) depending on parameters a and po- 



Values of 


1 + a > 0, 


1 + a > 0, 


1 + a < 0, 


1 + a < 0, 


H = a/a 


apQ > 


apo < 


apo < 


ctpo > 


H <-^ 


Contraction 




Contraction 




H = -^fpK 


Non-steady 


Expansion, 


Attractor 


Contraction 




equilibrium 


bounce and 




from the 




point 


contraction 




anti-Big Rip, 




Contraction, 




Expansion, 


bounce and 




bounce and 




bounce and 


expansion to 




expansion 




contraction 


the Big Rip 


H = ^ 


Attractor 




Non-steady 










equilibrium 










point 




H > ^ 


Expansion 




Expansion 





only restriction is the age of the universe which is hmited (e.g. by the oldest globular 
clusters) to tQ> 12- 10^ yr. We find the age of the universe in our model by integration 
of (23) from 2; = to 2; = oo for different values of a and po ^-nd for nowadays values 
of Ho = 65 km s"^ Mpc-\ Vtm = 0.3 and = 4.54 • IQ-^jh'^, where h is the Hq in 
the units 100 km s~^ Mpc~^. The connection of commonly used parameter w with our 
parameters a and po is given by relation w = a{l — Po/pde)- For example, if w = —0.82 
then the age restriction is important only at po/poE < —100. In the case of w < — 1 
the universe is older then in a pure A-term case {w = —1) and the age restriction is 
unimportant. 

In contrast, at Po/pde > 1 only a small range of parameters may correspond to the 
real universe. In addition to restriction on w and the universe age restriction one must 
require the universe to start from the Big Bang and not from the bounce at a finite 
scale factor. The bounce d = appears in the case of a > 1/3 and 0,^,0 < 0. Therefore, 
to avoid the bounce we must additionally require that a < 1/3. One more condition 
flafi + ^r,o > should bc satisfied in the boundary case of a = 1/3. Strictly speaking, 
the bounce at rcdshift 1 + Zb = {Qrfi/l^afilY^^^'^''^^ could not occur later then the end 
of inflation or reheating moment at the temperature Trh ~ 10^^ GeV corresponding 
to z-RH ~ 10^^. Therefore it is necessary to satisfy the condition Zb S> 2;rh- The last 
condition puts the restriction to the allowable values of a and fla,o- This restriction can 
be considered in the particular inflation model with taking into account the additional 
dynamical components (say inflaton). This consideration is out of scope of this paper. 




Figure 6. The evolution of da.rk energy witii an arbitrary equation of state p{p) is 
shown. The curves p = p{p) are the generalization of a linear equation of state, shown 
in Figures (2), (3), (4) and (5). Near each point of a sufficiently smooth curve p = p{p) 
one can represent the curve by its linear approximation thus reducing to the results of 
Section (2). 



4. Discussion and Conclusion 

In this paper we examined the dynamical evolution of the universe filled with a dark 
energy obeying the linear equation of state (3). It turns out that this simple linear 
model for the different choices of parameters a and po has a rich variety of cosmological 
dynamics. In dependence on signs of 1 + a and apo and the initial conditions for p and 
p there can be a set of distinctive types of the cosmological scenarios: Big Bang, Big 
Crunch, Big Rip, anti-Big Rip, solutions with the de Sitter attractor, bouncing solutions, 
and their various combinations. In the framework of the linear model (3) the analytical 
solutions of the dynamics of the universe were obtained. Using the phase plane analysis 
we gave the full classification of the solutions in dependence on the parameters a and 
Po- 

We distinguish four main types of the evolution of the universe filled with dark 
energy with a linear equation of state (3): 

(i) For parameters 1 + a > and apQ > the universe may contain cither non- 
phantom or phantom energy. For the non-phantom energy there arc two types of 
evolution: a) starting from the Big Bang the universe approaches to the de Sitter phase 
during an infinite time; b) reversed process with respect to the described above when 
the universe starts from a reversed in time the de Sitter phase and evolves to the Big 
Crunch. In the phantom case the universe starts from the reversed in time the de Sitter 
phase then it riches the bouncing point {p — 0, p ^ 0) and after bounce the universe 
approaches to the de Sitter phase. 

(ii) For parameters 1 -|- « > and apo < the universe may contain only non- 
phantom energy. In this case the only one type of evolution is possible: the universe 
expands starting from the Big Bang, reaches the bounce point {p = 0, p ^ 0) in finite 
time and then its expansion changes to the contraction, resulting in the Big Crunch. 
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(iii) For parameters 1 + a < and apQ < the universe may contain either non- 
phantom or phantom energy. In the case of phantom universe there arc two possible 
scenarios: a) starting from Oj = the universe expands in infinite time to the final 
Big Rip; b) the reversed in time process with respect to the described above when the 
universe starts from the anti-Big Rip and evolves to final state with a/ = 0. In the 
non-phantom case the universe starts from the reversed in time the de Sitter phase. 
Then it reaches the bouncing point (p = 0, p 7^ 0) and after the bounce the universe 
approaches the de Sitter phase. 

(iv) For parameters 1 -|- a < and ctpo > the universe may contain only phantom 
energy. The universe is born in the anti-Big Rip state and contracts reaching the 
bouncing point in a finite time. After the bouncing the universe begins to expand. In 
a finite time the universe comes to the Big Rip. 

The Table 1 summarizes the above results for the evolution of the universe filled 
with dark energy with a generalized equation of state (3). 

It should be stressed that the linearity of a considered dark energy equation of 
state is not crucial for the general properties of cosmological evolution. Instead of (3), 
one may consider any the rather smooth curves p — p{p) as shown in Fig. (6). It is 
clear that the general behavior of the evolution and properties of attractors and bounce 
points do not change in this case because any sufficiently smooth function p = p(p) can 
be linearized in the local vicinity of any point. Thus we can reduce the general problem 
for p = p{p) to the analysis of a linear cosmology considered in this paper. From the 
above it follows in particular that the universe filled with a dark energy with an equation 
of state = Pde(pde) always approaches the de Sitter attractor if additionally the 
physically reasonable conditions < ^Pde/c^Pde < 1 is satisfied. The first inequality in 
the last expression means that the considered dark energy is hydrodynamically stable. 
While the second inequality restricts the sound speed to the speed of the fight. A 
more detailed analysis of the cosmology with an arbitrary continuous equation of state 
p — p(p) will be presented elsewhere [20] . 

Our analysis is limited by the consideration of the cosmological model of the 
universe evolution filled only with dark energy. Taking into account the presence 
of the ordinary matter and radiation make the evolution of the universe much more 
complicated and is requied special consideration [20]. 
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